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Abstract 

We solve the spectrum of the closed Temperley-Lieb quantum spin chains us- 
ing the coordinate Bethe ansatz. These models are invariant under the quantum 
group U q [sl(2)]. 



1 Introduction 



Quantum group together with Temperley-Lieb algebra play a important role in the 
study of integrable spin chains. It may be interesting to study particular Hamiltonians 
associated with the Temperley-Lieb which are invariant to the quantum group. Taking 
into account usual toroidal boundary conditions, the Hamiltonians take the form 

N-l 

H= £ U n>n+1 + U m . (1.1) 

n=l 

where U ntn+ i operates in a direct product of two (2s + l)-dimensional complex spaces 
V 2s+l at positions n and n + 1. They are not invariant with respect to U q [sl(2)\ since 
ZTjvi 7^ U\n breaks translational invariance, reflecting the non-cocommutativity of the 
co-product. Indeed, we know from[|l], ||, ||] that very special boundary terms must be 
considered when we seek these quantum group invariant spin chains. In particular, 
one possibility to obtain a quantum group invariant Hamiltonian is to consider open 
boundary conditions, i.e., Uni = . For The XXZ-Hamiltonian with open boundary 
conditions one has to apply the Bethe ansatz techniques introduced by Sklyanin||] us- 
ing Cherednik's reflection matrices ||, ||. By this method the XXZ-Heisenberg model 0, 
the spl q (2, 1) invariant supersymmetric t-J model|8|, [| , the U g [sl(n)] invariant gener- 
alization of the XXZ-chainflOfl and the SU q (n\m) spin chains ||, |TTJ have been solved 



for open boundary conditions by this method. 

Recently, by means of a generalized algebraic nested Bethe ansatz, Karowski and 
Zapletal|n| presented a class of quantum group invariant n-state vertex models with 
periodic boundary conditions. Also an extension of this method to the case of graded 
vertex models was analyzed in [I3|, where a spl q (2\l) invariant susy t-J model with 



boundary conditions was presented. 

In fact, this type of models were first discussed by Martin [|TJJ from the representa- 
tions of the Hecke algebra. The study of closed quantum group invariant closed spin 
chains in the framework of the coordinate Bethe ansatz was presented by Grosse at al. 
for the SU q (2) case [I5|. In this context it would be interesting to discuss other quan- 



tum group invariant closed spin chains. Therefore, it is the purpose of this paper to 
present and solve, via coordinate Bethe ansatz [[nj a closed spin-s Hamiltonian, which 
in terms of the Temperley-Lieb operators can be written as 



N-l 



H=J2Un+Uo (1.2) 



n=l 



where 

U = GU^G- 1 , G = {Q-U X ){Q-U 2 )---{Q-U N „ 1 ) (1.3) 
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satisfying [H, G] = and additionally invariance with respect to the quantum algebra. 
The operator G shifts the U n by one unit GUnG^ 1 = U n+ i and maps U into U\, which 
manifest the translational invariance of Ti. 



2 The Temper ley- Lieb Hamiltonians 

In the basis where is diagonal with eigenvectors \s,n) ,\s — l,n) , \—s,n) and 
eigenvalues s, s — 1, —s, the Hamiltonian densities acting on two neighboring sites 
are given by 

(k,l\U\i,j) = e(i)e(k)q^5 i+jfl 5 k+lfi ^ ^ 

i,j,k,l — s, s — 1, • • • , — s + 1, — s 

where e(i) = (—1)* for s integer and e(i) = (— 1)* +1 / 2 for s semi-integer. Thus U n 
denotes the projection on states whose restriction to sites n and n + 1 has total spin 
zero. These Hamiltonians were derived Batchelor and KunibajlT]] from representations 
of the Temperley-Lieb algebras associated with quantum group U q [sl(2)}. The case 
s = 1/2 was investigated in reference Jl5 . 



In fact, U n obeys the Temperley-Lieb algebra] 18 



Ul = (Q + Q~ l )U n , Q + Q- 1 = [2s + 1], 
U n U n+1 U n = U n , [U w Ui] =0 for \n-l\>2 (2.2) 

and commutes with the quantum group U q [sl(2)]. The g-number notation is [x] q = 
(q x — q~ x )/(q — q" 1 )- This algebra appears in a large class of solvable models and 
is known to essentially govern their physical properties: Ti is an element of a set of 
infinity quantities conserved which are involutive provided that U n satisfies the defining 



relations (|2.2|) 



Having now built common ground for all closed Hamiltonian densities, whose salient 
feature is that they are spin-zero projectors, we may implement the steps of[[H|] , where 
the spectrum of the A-D Temperley-Lieb Hamiltonians with either periodic and free 
boundary conditions were solved , via a generalization of the coordinate Bethe ansatz. 



3 The coordinate Bethe ansatz 

Since these Hamiltonians commute with the total spin = Y.n=i ^n-> ^ ne eigenvalues 
of the operator r = sN — Sf, can be used to collect the eigenstates of TC in sectors, 
\l/ r . Due to this U(l) invariance, there always exists a reference state \l/o satisfying 
H^o = E ^o, with E = 0. We take ^ to be ^ = Un \s,n). This is the only 
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eigenstate in the sector r = 0. All other energies will be measured relative to this 
state. 

We will now start to diagonalize 7i in every sector. Nothing interesting happens 
in sectors with r < 2s. Since 7i is a sum of projectors on spin zero, these states are 
annihilated by TC. 

The first nontrivial sector r = 2s, the correspondent eigenspace is spanned by 
the states \n(— j,j)} = s s ■ ■ ■ s —jjs---s) , where n = 1, 2, N — 1 and j = 
—s,—s + l,...,s. We seek eigenstates of 7i which are linear combinations of these 
vectors. It is very convenient to consider the linear combination 

l«(n)>= E(-l) s+J ' Q s ' j H-j,j)). (3.1) 

j=-a 

which is a highest weight state, S + \ fi(n)) = 0, and eigenstate of U n 

U n \Q(n)) = (Q + Q- 1 )!^)), U n±1 \Q(n))=e s \Q(n±l)), 
U n \Q(n±l)) = e s \Q{n)), U n \n(m)} = for n^{m±l,m} (3.2) 

where e s = — 1 for s semi-integer and e s — 1 for s integer. In this basis, all spin-s 
Hamiltonians 7i can be treated in a similar way and it affords a considerable simplifi- 
cation in the diagonalization of 7i, when one compares with the computations in the 



usual spin basis JT9 



3.1 One-pseudoparticle eigenstates 

Let us consider one free pseudoparticle as a highest weight state which lies in the sector 
r = 2s 

N-l 

^2s = E A (n) |ft(n)) . (3.3) 

n=l 

Using the eigenvalue equation H^2s = E 2s ^2s, one can derive a complete set of equa- 
tions for the wavefunctions A(n). 

The action of the operator G = (Q — U\) ■ ■ ■ (Q — Un-i) on the states |fi(n)) can 
be computed using (|3.2|) . It is simple on the bulk and at the left boundary 



G \tt(n)) = -e s Q N ~ 2 |0(n + 1)) , 1 < n < N - 2 (3.4) 
but manifests its nonlocality at the right boundary 
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JV-1 



G \n(N - 1)> = e s Q N ~ 2 (-£sQ)- n \n(N - n)) (3.5) 

n=l 

Similarly, acting with the operator G^ 1 = (Q^ 1 — E/jv-i) ' ' ' (Q 1 ~ U\) , we get 

G' 1 |0(n)) = -e s Q~ N+2 \Q(n - 1)) , 2<n<N-l (3.6) 

2V-1 

G- 1 |0(1)> = e.g-" +a £ (-e s Q) n |0(n)> (3.7) 

71=1 

for the bulk including the right boundary and for the left boundary, respectively. 

From these results one can see that the action of Uq = GUn-iG^ 1 vanishes on the 
bulk 

W |O(n))=0 , 2<n<N-2 (3.8) 
and is nonlocal at the boundaries 

JV-l 

Uo = -e s £ (-tsQY |fi(n)) , Uo MN - 1)) = {-e s Q)- N U (3.9) 

n=l 

Next, the action of the operator W = Y,k=i Uk on the states \Q(n)) gives the 
following equations 

U\n{l)) = (Q + g- 1 )|fi(l)) + e s |fi(2)) 

W|fi(n)) = (g + g^ 1 ) |Q(n)) + e s - 1)) + e s |fi(n + 1)) 

for 2 < n < N — 2 

U\Q(N-1)) = (Q + Q- 1 )\Q(N-l)) + e s \fl(N-2)). (3.10) 

Before we substitute these results into the eigenvalue equation, we will define two 
new states 

e s |fi(0)) = W , e s \Q(N)) = U \Q(N - 1)) (3.11) 

to include the cases n = and n = N into the definition of ^2s, equation ( |3.3| ). Finally, 
the action of H = U + Uq on the states is 

w|fi(o)) = (g + g- 1 )|fi(o)) + (-e s g) iV e s |fi(iv-i)) + es |fi(i)) 

H\n(n)) = (Q + Q" 1 ) |fi(n)) + e a |fi(n - 1)) + e. |n(n + 1)) 
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for 1 < n < N — 2 

n\n{N-l)) = (Q + Q- 1 )\n(N-l)) + e s \n(N-2)) 
+ (-e s Q)- N e s \Q(0)) 

H\n(N)) = (Q + Q- 1 )|fi(iV))+e s |fi(iV-l)) 

+ (-e s Q)- N e s \Q(l)) (3.12) 

Substituting these results into the eigenvalue equation Tt^ 2s — E 2s and using the 
boundary conditions 

(-e s Q) N A(x) = A(N + x) (3.13) 
we get a complete set of eigenvalue equations for the wavefunctions 

E 2s A{n) = (Q + Q -1 )A(n) + e s A(n - 1) + t s A(n + 1) 

for 1 < n<N-l (3.14) 

The plane wave parametrization A(n) = £ n solves these eigenvalue equations and 
the boundary conditions (|3.13| ), provided that: 



E 2a = Q + Q- 1 + e s (Z + C 1 ) and £ N = (-e s Q) N (3.15) 
where £ = e td and 6 being the momentum. 

3.2 Two-pseudoparticle eigenstates 

Let us now consider the sector r = 2 (2s), where we have two interacting pseudopar- 
ticles. We seek the corresponding eigenfunction as products of single pseudoparticles 
eigenfunctions, i.e. 

^ 4s = E A{x 1 ,x 2 )\Q{x 1 ,x 2 )) (3.16) 

Xl+l<X2 

where 

\n( Xl ,x 2 )) = E E (-i) i+ V s -^>iH,0^2(-i,i)) (s.ir) 

i=~s j=—s 

To solve the eigenvalue equation H^4 S = E^^, we recall Q3.2|) to get the action 
of U and U on the states \Q(xi,x 2 )). We have to consider four cases: (i) When the 
two pseudoparticles are separated in the bulk, the action of U is 

U\n(x x ,x 2 )) = 2(Q + Q- 1 )\n(x 1 ,X2)) + e s \n(x 1 -l,x 2 ))+e s \n(x 1 + l,x 2 )) 

+e s \n( Xl ,x 2 - 1)) + e s \Q(xx, x 2 + 1)) (3.18) 
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i.e., for x\ > 2 and X\ + 3 < x 2 < N — 2; (ii) When the two pseudoparticles are 
separated but one of them or both are at the boundaries 

U\n(l,x 2 )) = 2(Q + Q- l )\n(l,x 2 )) + e s \n(2,x 2 ))+e s \n(l,x 2 -l)) 

+e s \Q(l,x 2 + l)) (3.19) 

U\Q(x u N-l)) = 2(Q + Q- 1 )\Q(x 1 ,N-l)) + e s \n(x 1 -l,N-l)) 

+e s \Sl( Xl + 1, N - 1)) + e s \Q(xi,N - 2)) (3.20) 

U N -l)) = 2{Q+Q- 1 ) AT - l))+e s |H(2, iV - l))+e s |0(1, N - 2)) (3.21) 

where 2 < x 1 < N — 4 and 4 < rr 2 < A 7 " — 2; (iii) When the two pseudoparticles are 
neighbors in the bulk 

U\Q(x,x + 2)) = 2(Q + Q- 1 )\n(x,x + 2)) + e s \n(x-l,x + 2)} + e s \n(x,x + 3)} 
+U x+1 \n(x,x + 2)) (3.22) 

for 2 < x < N — 4 and (iv) When the two pseudoparticles are neighbors and at the 
boundaries 

U 3)) = 2{Q + Q- 1 ) 3)) + e s 4)) + U 2 |0(1, 3)) (3.23) 

W|fi(iV-3,iV-l)) = 2(Q + g- 1 )|fi(A^-3,A^-l)) + e s |fi(A^-4,A^- 1)) 

+C^_ 2 |fi(JV-3,JV-l)) (3.24) 

Moreover, the action of Uq does not depend on the pseudoparticles are neither 
separated nor neighbors. It is vanishes in the bulk 

U \Q( Xl , x 2 )) = for x x ±\ and x 2 ^ N - 1, (3.25) 

and different of zero at the boundaries: 

X2-2 

U \n(l,x 2 )) = -e s ]T {-e s Q) k \Sl(k, x 2 )) - {-e^f^U^ \Q(x 2 -l,x 2 + 1)) 

k=i 

N-l 

-e s ]T (-e s Q) k - 2 Mx 2 ,k)) (3.26) 

k=x 2 +2 

Uo \n(x u N-l)) = (-e s Q)- N+2 Uo x 2 )) (3.27) 
where 2 < x 1 < N - 3 and 3 < x 2 < N - 2. 
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Before we substitute these expressions into the eigenvalue equation, we define new 
states in order to have consistency between bulk and boundaries terms 

U \n(l,x 2 )) = e s \Q(0,x 2 )), U \n(x 1 ,N-l)) = e a \n(x 1 ,N)) 
Uo\Sl(l,N-l)) = e.\n(p,N-l)) + e a \Sl(l,N)) 
U x+1 \n(x,x + 2)) = e s \tt(x,x + l)) + e s \tt(x + l,x + 2)) (3.28) 

Acting with TC on these new states, we get 

H\n(0,x 2 )) = 2(Q + Q- 1 )\Q{0,x 2 )) + e s \Q{0,x 2 -l)) + e s \Q{0,x 2 + l)) 

+e s x 2 )) + (-e s Q) N - 2 e s \Q(x 2 , N - 1)> (3.29) 

H\tt( Xl ,N)) = 2(Q + Q- 1 )\n(x 1 ,N)) + e s \n(x 1 -l,N))+e s \n(x 1 + l,N)) 

+e s \n( Xl , N-l)) + (-e s Q)- N+2 e s a*)) (3.30) 

H \n(x, x + 1) = (Q + Q- 1 ) \Q(x, x + 1) + e s \Q(x - 1, x + 1) + e s \Q(x, x + 2) (3.31) 

Substituting these results into the eigenvalue equation, we get the following equations 
for wavefunctions corresponding to the separated pseudoparticles. 

E As A( Xll x 2 ) = 2(Q + Q~ 1 )A(x 1 ,x 2 ) + e s A(xi - l,x 2 ) + e s A( Xl + l,x 2 ) 

+e s A(x u x 2 - 1) + e s A(x!,x 2 + 1) (3.32) 

i.e., for x±>l and X\ + 3 < x 2 < N — 1. The boundary conditions read now 

A(x 2 , N + Xl ) = {-e s Q) N - 2 A{x u x 2 ). (3.33) 

The parametrization for the wavefunctions 

A(x u x 2 ) = Aud 1 ^ 2 + A 21 £?Z? (3.34) 
solves the equation ( |3.32|) provided that 

E As = 2(Q + Q- 1 ) + e s (Ci + + 6 + (3-35) 

and the boundary conditions (|3.33|) provided that 

£ 2 = {-tsQ) -j— , §i = {-e s Q) — =>• £ = (-e s Q) ^ ; (3.36) 

^12 -^21 

where £ = £1^2 = e^ 6>1+6 ' 2 ' ) , 6*! + # 2 being the total momenta. 
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Now we include the new states ( |3.28| ) into the definition of \1/4 S in order to extend 

( PD to 

*4 S = E A(x u x 2 )\n(x u x 2 ). (3.37) 

Xl<X2 

Here we have used the same notation for separated and neighboring states. 
Substituting (|3.22|) and ( |3.31|) into the eigenvalue equation, we get 



E 4s A(x, x + 1) = (Q + Q^)A(x, x + 1) + e s A(x - 1, x + 1) + e s A(x, x + 2) (3.38) 

which gives us the phase shift produced by the interchange of the two interacting 
pseudoparticles 

A 21 1 + £ + e,(Q + 



A 12 l+£ + e,(Q + Q- 1 )£ 1 ' 
We thus arrive to the Bethe ansatz equations which fix the values of £1 and £2: 



(3.39) 



£2" = {-e s Q) N ~ 2 



1 + £ + e s (Q + Q- 1 )^ 



= (-^Q) 2(7V - 2) (3.40) 
3.3 General eigenstates 

The generalization to any r multiple of 2s is in principle straightforward. Since the 
Yang-Baxter equations are satisfied, there is only two-pseudoparticle scattering (if 
we use the S-matrix language). Therefore neighbor equations, where more the two 
pseudoparticles become neighbors, are not expected to give any new restrictions. For 
instance, in the sector r = 3 (2s) we have three interacting pseudoparticles with pa- 
rameters £i,£ 2 and £3. The corresponding wavefunctions 

A{x x , x 2 , x 3 ) = A l2 ^C 2 2 i? + A^Ci 1 ^ 3 ^ 2 + Am^fCs 3 + A 2al g'&& 

+^3i2£r£ 2 Cl £ 3 X2 + As^&Cz 1 (3.41) 

satisfy the boundary conditions 

A(x 2 , x 3 , N + xx) = (-e s Q) N ~ 4 A(x 1 ,x 2 ,x 3 ) 

which imply that 

^ = (-e s Q) N ~^ = (-e s Q) N -^, ^j^ = W (3.42) 
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These relations show us that the interchange of two pseudoparticles is independent of 
the position of the third particle. Thus in the sector r = p(2s), we expect that the 
p-pseudoparticle phase shift will be a sum of m two-pseudoparticle phase shifts and 
the energy is given by 

E P (2s) = £ {Q + Q- 1 + e s (£n + C 1 )} (3-43) 

n=l 

where 

& = (-e s Q) N ~ 2p+2 ft (- TZ7TZ € fn t 7^ I » a = W 
6 A / a 1 l + £,& + w + Q x )6 J 

= (-e,g) p(iV - 2p+2) (3.44) 

It is not all, in a sector r we may have p pseudoparticle and N s _i, N s _ 2 , N_ s+ i 
impurities of the type (s — 1), (s — 2), (— s + 1), respectively, such that 

N s „t + 2iV s _ 2 + • • • + (2s - 1)AL S+1 = r-2sp (3.45) 

We call impurity a state |a, n) flanked by at least two states \b, nil) such that a+b ^ 0. 
Since 7i is a sum of projectors on spin zero, these states are annihilated by Ti . In 
particular, the do not move under the action of 7i, which is the reason for their name. 
Nevertheless, a pseudoparticle can propagate past the isolated impurity, but in so doing 
causes a shift in its position by two lattice sites. Thus, for a sector r with I impurities 
with parameters £i, and p pseudoparticles with parameters £1+1, £i +p the energy 
is given by ( |3.43|) , and the Bethe equations do not depend on impurity type and are 
given by 



i+p 



gtkt -&= (-esQ) N ~ 2p+2 II - TT77 " nTn-i, ( 3 " 46 ) 



b=l+l 



l + Ub + esiQ + Q- 1 )^ 



with a = I + 1, 1 + 2, I + p , p > 1, and 

£ 2p (6 + i • • • M N -* = (~e s Q)^ N - 2 ^ (3.47) 
where £ = £i£ 2 • • • f • • • 

4 Conclusion 

We have shown that these closed Temperley-Lieb quantum invariant spin chains can 
be solved by the coordinate Bethe ansatz. A consequence of the nonlocal terms Uq is 
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the arising of boundary conditions depending on the quantum group parameter q via 
the relation Q + Q^ 1 = [2s + l] q and on the number p of pseudoparticles (which is 
equal to spin sector r , when s = 1/2). 

An interesting extension of this work would be the application of the methods 
here presented to solve new strongly correlated electronic systems associated with the 
Temperley-Lieb algebras [[2(], |21|]. This is presently under investigation 
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